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Abstract 

We present the results of an application of Bayesian inference in test- 
ing the relation between risk and return on the financial instruments. On 
the basis of the Intertemporal Capital Asset Pricing Model, proposed in 
|13| we built a general sampling distribution suitable in analysing this 
relationship. The most important feature of our assumptions is that the 
skewness of the conditional distribution of returns is used as an alterna- 
tive source of relation between risk and return. This general specification 
relates to Skewed GARCH-In-Mean model proposed in [14]. In order to 
make conditional distribution of financial returns skewed we considered a 
constructive approach based on the inverse probability integral transfor- 
mation presented in details in |15| . In particular, we apply hidden trun- 
cation mechanism, two equivalent approaches of the inverse scale factors, 
order statistics concept. Beta and Bernstein distribution transformations, 
and also the method recently proposed in [6|. Based on the daily excess 
returns on the Warsaw Stock Exchange Index we checked the empirical 
importance of the conditional skewness assumption on the relation be- 
tween risk and return on the Warsaw Stock Market. We present posterior 
probabilities of all competing specifications as well as the posterior anal- 
ysis of the positive sign of the tested relationship. 

PACS 89.65 Gh, 05.10 Gg 

1 Introduction 

The basics of the financial economics is constituted by the relationship between 
risk and return. Numerous papers have investigated this fundamental trade-off 

'Presented at the 3-rd Symposium on Socio- and Econophysics, FENS2007, Wroclaw 22—24 
November 2007. To be published in Acta Physica Polonica A 



1 



testing linear dependence of excess return on the level of risk, both measured by 
conditional mean and conditional standard deviation of the value of investor's 
wealth. According to [13], given risk aversion among investors, when investment 
opportunity set is constant, there is a positive relationship between expected 
excess return and risk. Hence, it is possible to express the risk in terms of the 
expected premium generated. 

Historically, authors have found mixed empirical evidence concerning the rela- 
tionship. In some cases a significant positive relationship can be found, in others 
it is insignificant and also some authors report it as being significantly negative. 
For instance, using monthly U.S. data [7] and also [3j found a predominantly pos- 
itive but insignificant relationship. Glosten, Jagannathan and Runkle reported 
in [8] a negative and significant relationship on the basis of Asymmetric-GARCH 
model, instead of commonly used GARCH-in-Mean framework; see j4j. Scruggs 
summarises the empirical evidence of considered relationship in [16] . 
The main goal of this paper is an application of Bayesian inference in testing the 
relation between risk and excess return of the financial time series. We revisited 
Intertemporal Capital Asset Pricing Model (ICAPM) in order to investigate the 
empirical importance of the skewness assumption of the conditional distribu- 
tion of excess returns. On the basis of the model, we built a general sampling 
distribution of the observables suitable in estimating risk premium. The most 
important feature of our model assumptions is that the possible skewness of 
conditional distribution of returns is used as an alternative source of relation 
between risk and return. Thus pure statistical feature is equipped with eco- 
nomic interpretation. Our general specification fully corresponds to suggestion, 
that systematic skewness is economically important and governs risk premium; 
see [To]. In order to make conditional distribution of financial returns skewed 
we considered a constructive approach based on the inverse probability integral 
transformation presented in details in [T^- Based on the daily excess returns of 
index of the Warsaw Stock Exchange we checked the total impact of conditional 
skewness assumption on the relation between return and risk on the Warsaw 
Stock Market. On the basis of the posterior probabilities and posterior odds 
ratios, we test formally the explanatory power of competing, conditionally fat 
tailed and asymmetric GARCH processes. 

2 Creating asymmetric distributions 

The unified representation of the univariate skewed distributions that we apply 
in this paper is based on the inverse probability integral transformation; for 
details see [T^. The family of random variables IP = {£s,£s : ^ M}, with 
representative density s{.\6,rip) is called the skewed version of the symmetric 
family I (of random variables with unimodal symmetric density f{-\0) and dis- 
tribution function F, such that the only one modal value is locaHsed at a; = 0) 
if s is given by the form: 

six\e,rjp)^f{x\e)-p{F{x\0)\r^p), xeR. (2.1) 



2 



The asymmetric distribution s{.\0,r]p), where 9 is inherited from the density / 
and r]p groups the skewness parameters, is obtained by applying the density 
pi-lrjp) as a weighting function. Within the general form l|2.ip several classes of 
distributions have been imposed on some specific famihes of symmetric random 
variables. A review of skewing mechanisms was presented in [15]. Here we 
apply l|2.ip in order to build the set of unnested specifications, which compete 
in explaining the possible relationship between risk and return. 

3 Basic model framework and competing skewed 
conditional distributions 

Let denote by Xj the value of a stock or a market index at time j. The excess 
return on xj , denoted by yj , is defined as the difference between the logarithmic 
daily return on Xj in percentage points {rj — 100 ln(a:j/a;j_i)) and the risk free 
short term interest rate (denoted by r*), namely yj = rj — r*. The volumi- 
nous literature focused on examination the relationship between risk and return 
bases on the Intertemporal CAPM, proposed by Merton in [13]. According to 
assumptions of Merton's theory there exists a set of distributions P(.|V'j-i), 
conditional with respect to the information set at time j (denoted by 
such, that: 

E{y,\^,^,)^a*D{y,\i,,^^), (3.1) 

where symbols E and D denote expectaction and standard deviation respec- 
tively. The coefficient a* > in p.ip measures the relative risk aversion of 
the representative agent. Under assumptions of the informational efficiency of 
the market, the information set at time j can be reduced to the history of the 
process of the excess return, namely V'i-i — {■■■! yj-2, Uj-i)- Consequently, an 
econometric model of the relationship between risk and return should explain 
the properties of the conditional (with respect to the past of the process yj) 
distribution of the excess return yj at time j. It is also of particular interest to 
find any linkage between expected excess return and the measure od dispersion 
of the distribution of conditional to V'j-i- Following [4], [7] and [T^ we 
consider for yj a simple GARCH-In-Mean process, defined as follows: 

yj^[a + E{zj)]^j+u,, j = l,2,..., (3.2) 

where Uj = [zj — E{zj)]^/hj, and Zj are independently and identically dis- 
tributed random variables with E{zj) < +oo. The scale factor hj is given by 
the GARCH(1,1) equation; see [2]: 

hj = ao + aiu'^_i + (3ihj-i, j = 1,2, 

The specific form of the conditional distribution of yj in (|3.2p is strictly depen- 
dent on the type of the distribution of Zj. Initially, in model denoted by Mq, we 
assumed for Zj the Student-t density with unknown degrees of freedom v > I, 
zero mode and unit inverse precision: 

Zj\Ma ^ iiSt{Q,l,v), v>l. 
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The density of the distribution of Zj in Mq is given as follows: 

2-1 -(i'+l)/2 



p(.iAfo) = /.(.io,i,.) = ^|§;±iii 



z 



Given model A/q, E{zj) — 0, Uj — Zj^fhj^ and hence p.2p reduces to simpler 
form ?/j = OL\f^j + Let denote by 6* = (a, ao, ai, the vector of all 
parameters in model Mq. Here, the conditional distribution of the error term 
is the Student-i distribution with degrees of freedom parameter u > 1, zero 
mode and inverse precision hj. Consequently, the following density represents 
conditional distribution of the excess return at time j: 

p(%|V,_i,0,A/o) = hj°-^, f tihj°-' {y, ~ a^j)\0.1, = 1,2,.... 

Given model Mq the expected excess return (conditional to the whole past V'i-i) 
is proportional to the square root of the inverse precision hj : 



E{yj\^,^i,e,Mo) = a^h,, j = l,2,.... (3.3) 

The parameter a G R captures the dependence between expected excess return 
and the level of risk, both measured by E{yj\ipj-i,9, Mq) and the scale param- 
eter ^yhj respectively. 

Now we want to construct a set of competing GARCH specicifations {Mi,i = 
1, . . . , fc} by introducing skewness into density of the conditional distribution of 
excess return, p{yj\ipj-i,0, Mq). The resulting asymmentric distributions are 
obtained by skewing the distribution of the random variable zj, according to 
method presented in the previous section. The asymmetric density of zj is of 
the general form related to the formula p.ip : 

piz\M,) = ft{z\0, 1, iy)p[Ft{z)\r,„M,],z e i?, i = 1, 2, . . . , fc, 

where p{.\'r]i,Mi) defines the skewing mechanism parameterised by vector rji 
and Ft{.) is the cdf of the standardised Student-i distribution. This leads to the 
general form of the conditional distribution of daily excess return yj in model 
A/,: 

p(y,|V',-i, 0,77., M,) = /iTO-5/t(z;|0, 1, iy)p[Ft{z*)\r,,,A'hlj = 1, . . . , 
where z* — hj'^'^{yj — Hj) and fij — [a + E{zj)]y/h'j. 

By imposing skewness, the expectation E{zj) is no longer equal to zero. Con- 
sequently, given Mi, the expectation of the excess return (conditional to V'j-i) 
is still proportional to ^/hj, but the coefficient of proportionality changes: 



E{yj I V-j- 1 , 0, 7]^ , M, ) = [a + £;(z, )] ^Jhj . 

Hence, the skewness of conditional distribution of yj is treated in Mi as another 
source of the tested relationship. 



4 



As the first specification, denoted by Mi, we consider GARCH model with 
skewed Student-t distribution obtained by the method proposed in [5], or equiv- 
alently in The model M2 is the result of skewing conditional distribution 
according to the hidden truncation idea; see In models M3 and M4 we apply 
Beta skewing mechanism; see citeJones2004. In model M5 we apply Bernstein 
density based skewing mechanism with m = 2 free parameters, while model 
Mg is built on the basis of the skewing construct defined in [6]. All compet- 
ing specifications, together with analytical forms of skewing mechanisms and 
model specific parameters are presented in Table 1. For some details concerning 
sampHng densities and prior specifications in each model see [TS] . 

4 Empirical results for WSE index 

In this part we present an empirical example of Bayesian comparison of all com- 
peting specifications. We also discuss the posterior analysis of the total impact 
of the conditional skewness assumption on the relationship between risk and 
return on the Warsaw Stock Exchange (WSE). Our dataset y was constructed 
on the basis of t=2144 observations of daily growth rates, rj, of the WSE in- 
dex (WIG) from 06.01.98 till 31.07.06. The risk free interest rate, r*, used to 
calculate excess return yj, was approximated by the WIBOR overnight inter- 
est rate (WIBORo/n instrument). Our empirical results remained practically 
unchanged for r* calculated on the basis of the middle and long term WIBOR 
Polish Zloty interest rates and also in the case r* ~ 0. 

Table 2 presents posterior probabilities P{Mi\y) calculated for each of competing 
models M^, ? = 0, 1, . . . , 6. The initial specification Mq, built on the basis of the 
conditional symmetric Student-i distribution, receives a little data support, as 
the posterior probability P{AIo\y) is slightly greater than 8%. All remaining pos- 
terior probability mass is attached to specifications which allow for conditional 
skewness. It is clear, that the modelled dataset of excess returns of WIG index 
do not support decisively superiority of any of the competing skewing mech- 
anism. The greatest value of P{Mi\y) receives conditionally skewed Student-t 
GARCH model generated by the Beta distribution transformation with two free 
parameters. In this case the value of posterior probability is equal about 40%. 
The dataset also support conditionally skewed Student-i GARCH model with 
hidden truncation mechanism (M2) and Beta distribution transformation with 
one free parameter (M3). Those three models cumulate more than 90% of the 
posterior probability mass, making all remained conditionally skewed specifi- 
cations improbable in the view of the data. Thus, inverse scale factors, the 
Bernstein density transformation and construct proposed in namely models 
Ml, M5 and Mq, lead to very doubtful explanatory power. Those specifications 
are strongly rejected by the data, as the values of posterior probabilities are 
much smaller than posterior probability of symmetric GARCH model. 
In Table 2 we also compare the total impact of the conditional skewness effect on 
the tested relation between risk and return. According to our assumptions, the 
conditional expectation of the excess return is proportional to the square root 
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of the inverse precision hj . Since we parameterize the market risk by a more 
general dispersion measure than standard deviation we report the information 
about the relative risk aversion by the posterior probability of the positive sign 
of the function a + E{zj). According to l|3.2p it enables to test the positive sign 
of the relative risk aversion coefHcient. Initially we checked the strength of the 
relation in model Mq, which does not allow for conditional skewness. Given Mq 
E{zj) = and the whole information about relative risk aversion is reflected 
in parameter a; see l|3.3p . Just Hke many researchers we obtained positive, but 
rather weak, relation between expected excess return and risk, given model with 
symmetric conditional distribution. The posterior probability P{a > 0\Mo,y), 
equal about 0.92, leaves considerable level of uncertainty about true strength of 
the tested relation. Consequently, model Mq does not conflrm our hypothesis 
strongly. Imposing unreasonable (in the view of the data) skewness into condi- 
tional distribution of excess returns also may not strengthen our inference. In 
case of models with weak data support (Mi and Mg) the assumption of asym- 
metry of the density p{yj\ipj-ij9,r]i, Mi) does not improve posterior inference 
about the sign of a -I- E{zj). In case of Mi and Mg posterior probability of pos- 
itive relationship is very close to the value generated within Mq. Only in case 
of the skewing mechanisms with the greatest data support, namely Beta trans- 
formation with two parameters and hidden truncation, the WIG excess returns 
yield decisive support of the positive sign of the relative risk aversion coefHcient. 
In case of model M3, the posterior probability of positive sign of a -I- E{zj) is 
greater than 0.99, leaving no doubt about the signiflcance of the relationship 
between risk and return postulated by Merton in [13]. Hence, it was possible 
to confirm positive sing of a -I- E{zj) strongly only by imposing specific skewing 
mechanism into conditional distribution of excess returns. Beta distribution 
transformation with two free parameters was able to detect additional source of 
information about risk premium in the WIG dataset. Also, hidden truncation 
mechanism and Bernstein density transformation strongly confirm positive sing 
of the risk aversion coefficient, as posterior probability P{a -\- E{zj) > 0|Mi,y) 
is greater than 0.98, for i = 2 and 5. 



5 Concluding remarks 

We checked the impact of the conditional skewness assumption on the strength 
of the relationship between risk and expected return. On the basis of the In- 
tertemporal CAPM model, proposed in [13], we built a GARCH-In-Mean type 
sampling distribution suitable in modelling such relationship. Our approach, 
which fully relates to the model proposed in [TJ, treats the skewness of the 
conditional distribution of excess returns as an alternative source of informa- 
tion about risk aversion. Based on the daily excess returns of the Warsaw Stock 
Exchange index we checked the empirical importance of the conditional skew- 
ness assumption on the relation between risk and return. Posterior inference 
about skewing mechanisms showed positive and decisively significant value of 
the coefHcient of the relative risk aversion once a specific skewing mechanism 
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was imposed in conditional Student-i distribution. The greatest data support, 
and also very strong support of the relation postulated by Merton in [13], re- 
ceived skewness generated by Beta distribution transformation with two free 
parameters. 
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Table 1: The conditional skewing mechanisms p{.\'r]i), defined for u £ (0,1), 
skewness parameters rji and conditional symmetry restrictions in all competing 
specifications Mj, i — 1,2, 3, 4, 5, 6. 



Ms 

Bernstein density (2 parameters) 

(UJI, UJ2) £ (0, 1)^, UI3 = 1 — LUl — UJ2 

p{u\uJi,W2) = J2'i-i ^jBe(u\ j, 4 - j) 
1 

symmetry: uji — uJ2 — — 


Ml 

Inverse scale factors; see \E\ and [9] 
parameterisation as in [5] 
2 

71 > 0, C = r- 

71 + 7i 

/l(tl) = /(0,0.5)("), hiu) = /[0.5.1)(") 

, , , ^f{jiF-^(u))hiu) + fij-^F-Hu))^^ 

''^"'-^^^^'^ /(F-W) 
symmetry: 71 — 1 


M2 


M3 


Hidden truncation; see [T] 


Beta one parameter; see ffl| 


72 £ R 


73 > 


p(u|72) = 2F(72F-H")) 


Piu\l3) ^ -Be(u|73, 7~^) 


symmetry; 72 — 


symmetry: 73 — 1 


Ms 

The construct proposed in [6] 

74 £ R 

p(u|74) = 1 + i(74)[3(M|74) - 1] 

symmetry; 74 — 


M4 

Beta two parameters; see and |l2] 
a> 0,b> 
p{u\a,b) — Be{u\a, b) 
symmetry: a — b — 1 



Table 2: Decimal logarithms of the marginal data density values, posterior 
probabilities of all competing specifications Mi and posterior probabilities of 



the positive sign of the relative risk aversion coefiicient a + E{zj). 





Ml 


Ma 


Ma 


M4 


Ms 


Me 


Mo 


log p{y\Mi) 


-1559.45 


-1558.50 


-1558.78 


-1558.41 


-1560.82 


-1560.10 


-1559.06 


P(Mi\y),i = 0, . . . , 6 
P(Mi|y),i = 1, . . . , 6 


0.0.353 
0.0387 


0.3152 
0.3452 


0.1654 
0.1811 


0.3878 
0.4246 


0.0015 
0.0017 


0.0079 
0.0087 


0.0868 


P(a + E{zj) > 0\Mi,y) 


0.9102 


0.9894 


0.9528 


0.9972 


0.9893 


0.9230 


0.9201 
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